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We characterize good clocks, which are naturally subject to fluctuations, in statistical terms. We 
also obtain the master equation that governs the evolution of quantum systems according to these 
clocks and find its general solution. This master equation is diffusive and produces loss of coherence. 
Moreover, real clocks can be described in terms of effective interactions that are nonlocal in time. 
Alternatively, they can be modeled by an effective thermal bath coupled to the system. 



PACS: 03.65.Bz, 05.40.+j, 42.50.Lc 
I. INTRODUCTION 

In quantum mechanics, spatial positions are described 
by quantum observables. However, the situation with re- 
spect to the time parameter is rather more involved, as 
has been known since Pauli pointed out that there is no 
self-adjoint operator canonically conjugate to the total 
energy, if the spectrum of the latter is bounded from be- 
low Q. Some attempts have been made to circumvent 
this problem by making use of the concept of time of ar- 
rival Recent developments, however, have mapped 
the problem of constructing a good time of arrival opera- 
tor into the problem of constructing a good time operator 
|| , thus showing that the special status of time is deeply 
embedded into the structure of quantum mechanics. We 
are then bound to use real physical clocks and rely on 
their readouts when measuring the evolution of a quan- 
tum system. 

Any real clock is inevitably subject to quantum fluc- 
tuations, which introduce uncertainties in the equations 
of motion. For instance, it has been shown that the fi- 
nite mass and size of the clock impose limitations in the 
measurement of spacetime distances in the framework 
of general relativity [Q. Some considerations have also 
been made on the role of quantum clocks in the context 
of quantum cosmology [j5|,^| . Simple models for quantum 
clocks have been proposed, and the quantum evolution of 
a system according to a quantum clock suitably coupled 
with it has also been studied || [|. The general con- 
clusion is that the system becomes more and more per- 
turbed as the resolution of the clock is improved. Even 
more, quantum gravity may well imply the existence of 
an absolute limit, the Planck scale, to the accuracy of 
spacetime-distance measurements and, in particular, to 
clock synchronization (for a review, see, e.g., Rcf. |lC|]), 
with possible effects in the low-energy regime Q . 

It follows that any quantum clock that we could pos- 
sibly build would lead to uncertainties and errors. These 
quantum errors, however, are not the only source of ran- 
domness in the measure of time. Real clocks are also 
subject to classical imperfections, small errors, that can 
only be dealt with statistically. For instance, an un- 
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avoidable classical source of stochasticity is temperature, 
which will introduce thermal fluctuations in the behav- 
ior of real clocks. Although this is not necessarily the 
most important source of errors in modern-day atomic 
clocks, it is nonetheless always present to some extent. 
In other words, the third law of thermodynamics forbids 
the existence of ideal clocks. 

We will study, within the context of the standard quan- 
tum theory, the evolution of an arbitrary system accord- 
ing to a real clock. Unlike other works we are 
not concerned about the quantum dynamical behavior 
of the clock but only with the readings of time that it 
provides. As stated before, these readings will undergo 
errors, which will be described by a stochastic process. 
It is in this sense that we regard these clocks as classi- 
cal, although this does not preclude a possible quantum 
mechanical origin for the stochastic fluctuations. Notice 
also that we will mostly consider clocks which are decou- 
pled from the system under study. In what follows we 
shall not delve further into the source of stochasticity, 
but assume a phenomenological description of it. 

Let us be more specific about the meaning of this ran- 
domness in the readings of the clock. Imagine a large 
ensemble of identical systems, prepare one of them in a 
given initial state at initial clock time t — 0, and then 
measure the state of that system at clock time t. If we 
repeat this procedure for all the systems in the ensem- 
ble, the result will be a probability distribution for the 
possible outcomes, its dispersion partially being a conse- 
quence of the lack of knowledge of the precise ideal time 
that has elapsed. Therefore, the evolution according to 
the readouts of the real clock is nonunitary. In other 
words, the use of real clocks induces loss of coherence in 
most physical quantum states, as we will explicitly show. 

To attain this objective we will start by analyzing in 
some detail the characteristics that should be expected 
from a good real clock. This is expounded in the next 
section. The implications of the randomness in the read- 
ings of a real clock f or t he quantum evolution of a system 
is explored in Sec. 
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lution according to real clocks from an effective point of 
view in terms of nonlocal interactions in time and also 
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in terms of a thermal bath coupled to the system under 
study. Dccohcrence owing to clock errors is investigated 
m Sec. ^ Finally, we summarize and conclude in the 
last section with some general remarks. 

II. GOOD CLOCKS 

Let us consider the phase space of a classical physical 
system, divided in a set of ordered cells, such that the 
expected evolution of the system moves the state from 
one cell fc to the next fc + 1 in a time e approximately. An 
observer keeps tally of which is the cell of largest k that 
has been visited, and thus the passage of time is recorded. 
Let Sk denote the ideal time at which the readout of the 
clock is ke, i.e., Sk is the ideal time of first arrival of 
the system at the fcth cell. For an ideal clock, the time 
s k and the clock time ke would coincide. However, for 
a real clock, these two quantities will differ by an error 
Afe = s k — ke, where it should be noticed that the index 
k pertains to the readout ke of the clock, i.e., to the 
fcth tick, and not to a preset ideal time. Given any real 
discrete clock, its characteristics will be encoded in the 
probability distribution for the stochastic sequence fl2l 
of clock errors, P({Afc}), which must satisfy appropriate 
conditions, so that it describes a good clock. 

A first property is that Galilean causality should be 
preserved, i.e., that causally related events should always 
be properly ordered in clock time as well, which implies 
that Sk+i > Sk for every k. In terms of the discrete 
derivative a k = (A k +i — Ak)/e of the stochastic sequence 
{Afc}, we can state this condition as requiring that, for 
any realization of the stochastic sequence, a k > — 1. 

A second condition that we would require good clocks 
to fulfill is that the expectation value of relative errors, 
determined by the stochastic sequence be zero, i.e., 

(ctk) = for all fc. If this were not the case, the clock 
would either systematically go fast or slow down, and a 
redefinition through this systematic drift would provide 
us with a well-centered clock. Consequently, the expecta- 
tion value for the absolute errors {A^} will be constant. 
Furthermore, since fc ~ will be the time at which the 
systems whose evolution we are studying are prepared, 
Ao will not be stochastic and, without loss of generality, 
will be set to zero by a simple translation of the origin of 
time, so that (A*,) = 0. 

Another expectation that we would have for a good 
clock is that it should always behave in the same way (in 
a statistical sense) . The difference between the ideal time 
for the (fc + l)th tick and the ideal time for the fcth tick 
must be always the same, statistically speaking, even if 
the actual errors A k +i and A k are large by some statisti- 
cal fluctuation. Therefore, we can say that the clock be- 
haves consistently in time as a good one if those relative 
errors {a k } are statistically stationary, i.e., the probabil- 
ity distribution V({ak}) for the sequence of relative er- 
rors {a.k\ (which can be obtained from P({Ak}), and vice 



versa) must not be affected by global shifts fc — > fc + fco of 
the readout of the clock. Note that the stochastic process 
A*, need not be stationary, despite the stationarity of the 
process a k - 

It is also intuitively obvious that the one-point proba- 
bility distribution function for the variables ak should be 
highly concentrated around the zero mean, if the clock 
is to behave nicely. Even more, it is to be expected for 
clocks with small errors that all the higher-order cumu- 
lants be much smaller than the correlation, which, in 
turn, should also be bounded by a small number, i.e., 

(a*) = , (a k a k -j) = Cj < cq < 1 , 

where Cj = C-j. The correlation for the sequence of 
absolute errors {A k } can then be easily obtained and 
has the form (A k Ai) = e 2 £\ =1 J2j=i c i-j- 

The correlation time $ for the stochastic sequence {a k } 

is given by ■& = £j2^-oo c i/(^ c o)- We wm introduce 
a new parameter n with dimensions of time, defined as 
ft 2 = coi? 2 . This comes about because, when the errors 
of the clock have a thermal origin, n 2 is proportional to 
the temperature, and independent of ?9. In general, the 
good-clock conditions imply k <C As we shall see, f? 
cannot be arbitrarily large, and, therefore, the ideal clock 
limit is given by n — » 0. 

Until now we have discussed general properties that 
a good clock must fulfill, regardless of the physical sys- 
tem under study. In addition to these properties, a good 
clock must have enough precision in order to measure the 
evolution of the specific system, which imposes further 
restrictions on the clock. On the one hand, the char- 
acteristic evolution time £ of the system must be much 
larger than the correlation time $ of the clock. On the 
other hand, the leading term in the asymptotic expansion 
of the variance (A 2 ) for large fc is of the form ^{ke/d) 
which means that, after a certain period of time, the ab- 
solute errors can be too large. The maximum admissible 
standard deviation in Ak must be at most of the same 
order as (. Then the period of applicability of the clock 
to the system under study, i.e., the period of clock time 
during which the errors of the clock are smaller than the 
characteristic evolution time of the system is approxi- 
mately equal to C 2 #/k 2 . For a good clock, k « i) « (, 
as we have seen, so that the period of applicability is 
much larger than the characteristic evolution time C- 

Even though, so far, we have only spoken of good 
discrete clocks, by analogy, we will consider continuous 
stochastic processes a(t), with corresponding probability 
functional V[a(t)]. The conditions previously stated for 
the discrete sequence {ak} admit a straightforward gen- 
eralization to the continuous case. In what follows, we 
shall use the formulation in the continuum. 
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III. MASTER EQUATION 

We shall now obtain the evolution equation for the 
density matrix of an arbitrary quantum system in terms 
of the clock time t. Let p s (s) be the density matrix for 
a quantum system whose unitary evolution in the ideal 
Schrodinger time s = t + A(t) is provided by the time- 
independent Hamiltonian H. 

For any given realization of the stochastic process that 
characterizes a good clock, and using the chain rule, we 
can write the von Neumann evolution equation in terms 
of the clock time t as 

d tPs {t + A(i)) - -t(l + a(t)) [H, Ps {t + A(t))] , 

where % has been set to 1. 

Let us now transform to the interaction picture in 
which the density matrix has the form 



pl(t + A(t))=e* m p s (t + A(t))e 



iHt 



Notice that the interaction term a(t)H has the same form 
in both pictures because it is proportional to the free 
Hamiltonian H. 

Integrating the resulting equation between and t, and 
reintroducing the result for p\, we obtain the following 
integro-differential equation: 

d t pl(t + A(t)) = -ia(t)[H,pl(0)} 



d£a{t)a{t') [H, [H, p\{t' + A(i'))] ] 



In order to find the evolution equation in the time t, we 
have to average this equation over all possible realizations 
a(t) of the stochastic process with the functional weight 
■p[a(£)]. The average of the density matrix p\(t + A(i)) 
will be denoted by p l {t) and can be regarded as the den- 
sity matrix of the system at clock time t. 

At the real time t = 0, we impose the initial condi- 
tion p : (0) = pl(0) = Pq. Additionally, for a good clock, 
(a(t)} — 0, as already discussed, and, as a consequence, 
the average of the linear term in a(t) vanishes. Further- 
more, the clock time derivative dt and the average over 
a(t') commute because V[a(t')] is stationary. Finally, the 
density matrix p\(t' +A(t')) can be expanded in powers of 
A(t'). Then the average of the integro-differential equa- 
tion for the density matrix p 1 yields 



P\t) 



[ dTc(r)[H, [H,p\t-T)]\ 
Jo 



0«a 3 )) 



where the overdot denotes derivative with respect to the 
clock time t. We have also performed a change of the 
integration variable from t' to r = t — t' and have in- 
troduced the correlation function c(r) for the stochastic 
process a(t). 

For a good clock, the higher order terms in a can be 
seen to be much smaller than the c(r) term by a fac- 
tor (k/C) 2 *C 1, provided that the system evolves for 



a time smaller than the period of applicability of the 
clock. Since £ 3> the system does not evolve signifi- 
cantly within a correlation time, and we can substitute 
p l (t — t) by p l (t). This is the so-called Markov approxi- 
mation. The process A(t) will not be Markovian in gen- 
eral and there is no reason for requiring that the process 
a(t) has this property either. However, and even though 
the Markov approximation refers to the system and not 
to the clock itself, it renders the possible non-Markovian 
character of the clock irrelevant. Furthermore, for evolu- 
tion times t much larger than the correlation time i?, we 
can take the upper integration limit to infinity. 

The resulting master equation, once we go back to the 
Schrodinger picture, can be written as 

p(t) = -i[H,p{t)] - (k 2 /0)[H, [H,p(t)]] . 

Notice that, in the ideal clock limit, n — * 0, the uni- 
tary von Neumann equation is recovered. We should also 
point out that this master equation is not a truncation 
of the BBGKY hierarchy fl|, and that irreversibility ap- 
pears because the errors of the clock cannot be eliminated 
once we have started using it. 

Under the good-clock conditions, k « i), we can ap- 
proximate "P by a stationary Gaussian probabil- 
ity functional with zero mean and correlation given by 
the correlation c(t) of V[a(t)}. Although this Gaus- 
sian approximation assigns a nonvanishing probability 
to a(t) < —1, this probability will be negligibly small 
since, for good clocks, c{t) <C 1. Thus the Gaussian ap- 
proximation to good clocks fulfills the Galilean causality 
condition for all practical purposes. 

In the Gaussian approximation, there is essentially 
only one good clock for which a(t) is Markovian, the 
Ornstein-Uhlenbeck process 12 1. In this case, the cor- 



relation function for a(t) in the stationary regime is 
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Since the possible non-Markovian 



C{T) 

character of the clock does not influence the time evolu- 
tion of the system (provided that the condition ( -C $ 
is satisfied, as happens for good clocks), the Ornstein- 
Uhlenbeck clock is generic in what concerns the evolution 
of quantum systems according to real clocks. 



IV. EFFECTIVE DESCRIPTIONS 

The master equation corresponds to the evolution of a 
system with a free Hamiltonian H coupled with a clas- 
sical noise source a(t), with a probability functional dis- 
tribution 'P[Q;(i)], via the interaction Hamiltonian a(t)H. 
The path integrals for this system then follow the pattern 
111 



VaV[a{t')] / VqVpe 



i[f dt(pq-H)-J dta(t)H(t)\ 



In the good-clock approximation, only the two-point cor- 
relation function c(r) is relevant, so that we can write the 
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probability functional as a Gaussian distribution. The 
integration over a(t) is then easily performed to yield 

J VqDpe I dt (?*i- ff M / dtdt'c(t-t')H(t)H(t') _ 

Therefore, we see that the effect of using real clocks for 
studying the evolution of a quantum system is the ap- 
pearance of an effective interaction term in the action 
integral which is bilocal in time. This can be understood 
as the first term in a multilocal expansion, which cor- 
responds to the weak-field expansion of the probability 
functional around the Gaussian term. 

This nonlocality in time admits a simple interpreta- 
tion: correlations between relative errors at different in- 
stants of clock time can be understood as correlations be- 
tween clock-time flows at those clock instants. The clock- 
time flow of the system is governed by the Hamiltonian 
and, therefore, the correlation of relative errors induces 
an effective interaction term, generically multilocal, that 
relates the Hamiltonians at different clock instants. 

The expression for the influence functional is, in the 
weak noise limit, completely analogous to the expres- 
sions above. Since the noise source is classical, we see 
that there is no dissipative term there, nor in the mas- 
ter equation [JL4J. Moreover, as the interaction term is 
proportional to H, there is no response of the system 
to the outside noise, which means that the associated 
impedance is infinite Jl5| , ^6| . 

From a different point of view, the clock can be ef- 
fectively modeled by a thermal bath, with a tempera- 
ture T to be determined, coupled to the system. Let 
H + Hi nt + -Hb be the total Hamiltonian, where H is the 
free Hamiltonian of the system and H D is the Hamilto- 
nian of a bath that will be represented by a collection of 
harmonic oscillators The interaction Hamiltonian 

will be of the form H lnt — £if , where the noise operator 
£ is given by 

V ^ Jo 

In this expression, a and are, respectively, the annihi- 
lation and creation operators associated with the bath, 
and x( w ) is a re &l function, to be determined, that repre- 
sents the coupling between the system and the bath for 
each frequency u>. 

Identifying, in the classical noise limit, the classical 
correlation function of the bath with c(r), the suitable 
coupling between the system and the bath is given by 
the spectral density of fluctuations of the clock: 

poo 

k B T h x(^) 2 = / drc(r) cos(wt) , 
Jo 

where k B is Boltzmann's constant. With this choice, the 
master equation for evolution according to real clocks is 
identical to the master equation for the system obtained 
by tracing over the effective bath. 



To go beyond the classical noise limit requires the in- 
troduction of the usual quadratic dissipation term in the 
influence functional p[ . However, the peculiar coupling 
£H implies that this term does not produce dissipation 
in the equations of motion: the fluctuation-dissipation 
theorem, which reflects the microscopic structure of the 
bath, is thus fulfilled, but there is no dissipation. 

V. DECOHERENCE 

The master equation contains a diffusion term and will 
therefore lead to a loss of coherence 0. However, this 
loss depends on the initial state. In other words, there ex- 
ists a pointer basis |I8j , so that any density matrix which 
is diagonal in this specific basis will not be affected by 
the diffusion term, while any other will approach a diago- 
nal density matrix. The stochastic perturbation a(t)H is 
obviously diagonal in the basis of eigenstates {\n}} of the 
Hamiltonian, which is therefore the pointer basis: the 
interaction term cannot induce any transition between 
different energy levels u> n . 

The components of the density matrix in this basis 
are p nm = (n\p\m). The master equation can be solved 
exactly, its general solution being 

U<) = /w(0)e— t e-('"»» , 

where u) nm — lu„ — uj m . The smallest energy difference 
to provides the inverse of the characteristic time for the 
evolution of the system, £ = 1/ui. The smallest decay 
constant is ui 2 k 2 /'0, equal to the inverse of the period 
of applicability of the clock. By the end of this period, 
the density matrix will have been reduced to the diago- 
nal terms and a much diminished remnant of those off- 
diagonal terms with slow evolution. In any case, the von 
Neumann entropy grows if the density matrix is not ini- 
tially diagonal in the energy basis. 

The effect of decoherence due to errors of real clocks 
does not only turn up in the quantum context. Consider 
for instance a classical particle with a definite energy E 
moving under a time- independent Hamiltonian H. Be- 
cause of the errors of the clock, we cannot be positive 
about the location of the particle in its trajectory on 
phase space at our clock time t. Therefore we have an 
increasing spread in the coordinate and conjugate mo- 
mentum over the trajectory. For a generic system, this 
effect is codified in the classical master equation 

q={H,q} + (k 2 /#){H,{H,q}}, 

where g(t) is the probability distribution on phase space 
in clock time. This classical master equation can be de- 
rived in a manner completely analogous to the quantum 
one. 

For simplicity, let us study the particular example of 
one-dimensional Hamiltonian motion with closed orbits, 
with H = luJ, ip being the angle variable with period 2it 
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conjugate to the action variable J, and w a constant fre- 
quency characteristic of the system. The classical master 
equation for the probability density g(tp, J; t) reads 

d t Q = ud vQ +{u J 2 K 2 /V)dlQ. 

This diffusion equation can be exactly solved by sepa- 
ration of variables. The slowest decaying mode has, as 
before, a decay constant uj 2 k 2 /§. 

In the case of one particle that is released with en- 
ergy E and initial angle ip$, the probability distribution 
spreads out over the corresponding connected component 
of the energy shell, and tends to 8{J — E/u>)/2w as clock 
time grows. As we can see, the information about the tp 
variable is washed out by the errors in our clock, which 
is precisely the information that is not available in the 
quantum case: if J is completely known for a given quan- 
tum state, the indeterminacy in its conjugate variable 
will be infinite, the situation towards which the classical 
decoherence process tends. 

Finally, it should be observed that the mechanism 
of decoherence is neither tracing over degrees of free- 
dom, nor coarse graining, nor dephasing |l7| , |l9| ]. Even 
though there is no integration over time introduced here 
by fiat, as happens in dephasing in quantum mechanics, 
the spread in time due to the errors of the clock has a 
similar effect, and produces decoherence. 



VI. CONCLUSIONS 

In our study of the evolution of quantum systems ac- 
cording to real clocks, which are necessarily subject to 
errors, we have first established a stochastic character- 
ization of good real clocks. Using this description, we 
have derived a master equation for the quantum evolu- 
tion in real-clock time and we have also found its general 
solution in the basis of energy eigenstates. The stochastic 
features of good real clocks and their effects on the quan- 
tum evolution can be equivalently described by means of 
interactions which are nonlocal in time. They can also 
be effectively modeled by a quantum thermal bath. The 
master equation exhibits a diffusion term which is re- 
sponsible for the loss of coherence of most initial states. 
Finally, we have analyzed the evolution of classical sys- 
tems according to real clocks and reached analogous con- 
clusions. 

The third law of thermodynamics and the quantum 
fluctuations prevent real clocks from being perfectly 
accurate. This suggests that, strictly speaking, the 
Schrodinger unitary evolution equation is just an exce- 
lent approximation valid for sufficiently short periods of 
time and that should be substituted, along the lines pro- 
posed in this paper, by a diffusive master equation in 
more general situations. This adds a random aspect to 
the evolution of quantum systems. Indeed, coherence is 
progressively lost until we reach the period of applica- 
bility of the clock and, after that time, unpredictability 



sets in, as we have seen. Even perfectly isolated systems 
will suffer loss of coherence because of the fluctuations of 
the real-clock and will appear as effectively coupled to a 
reservoir. 
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